A Tale of Two Theories: Quantum Griffiths Effects in Metallic Systems 
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We show that two apparently contradictory theories on the existence of Griffiths- McCoy singulari- 
ties in magnetic metallic systems JJ0 are in fact mathematically equivalent. We discuss the generic 
phase diagram of the problem and show that there is a non-universal crossover temperature range 
T* < T < luo where power law behavior (Griffiths-McCoy behavior) is expect. For T < T* power 
law behavior ceases to exist due to the destruction of quantum effects generated by the dissipation 
in the metallic environment. We show that T* is an analogue of the Kondo temperature and is 
controlled by non-universal couplings. 
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The problem of non-Fermi liquid behavior (NFL) in U 
and Ce intermetallics continues to attract a lot of atten- 
tion due to the breakdown of Landau's Fermi liquid the- 
ory in metallic alloys NFL is often characterized by 
power law or logarithmic temperature behavior in phys- 
ical quantities such as the magnetic susceptibility, x(T) 
and specific heat, Cy(T). Many NFL materials have as 
common features the closeness to a magnetic phase tran- 
sition and disorder generated by the alloying |J]. It has 
been proposed that NFL behavior can be associated in 
some materials with quantum Griffiths-McCoy singular- 
ities close to a quantum critical point (QCP) 0- These 
singularities are related with the tunneling, at low tem- 
peratures, of magnetic clusters with N spins generated by 
the percolating nature of the magnetic phase transition. 
As a result, the physical properties acquire non-universal 
power law behavior, X {T) oc C V {T)/T oc T~ 1+x , with 
the exponent A < 1 dependent on the distance from the 
QCP. While in its first version [5( the theory did not con- 
sider properly the dissipation coming from the electronic 
degrees of freedom, it was extended to include dissipa- 
tion [jj with the final conclusion that power law behav- 
ior disappears below a crossover temperature T* due to 
the freezing of the magnetic clusters (that is, the quan- 
tum Griffiths singularities are suppressed close enough 
to the QCP). Instead of power law, it was predicted 
that for T < T* one should have x(T) oc l/(Tln(l/T)) 
and C V (T)/T oc l/(Tln 2 (l/T)) (sec figQJ. These re- 
sults seem to explain a stronger than power law diver- 
gence and also the discrepancy between the specific heat 
and susceptibility exponents at low temperatures in cer- 
tain materials |fj. Nevertheless, for T* < T < luq it 
was found that x(T) oc T~ 1+A , that is, quantum Grif- 
fiths behavior. Besides the temperature dependence the 
theory has also predicted a magnetic field, H, depen- 
dence of the magnetization, M(H, T), and Cy(H, T) with 
M(H,T) oc HT- 1+X and C V (H,T)/T oc T~ 1+x for 
T > (H,T*) and M(H, T) oc H x and C V (H,T)/T oc 
T -3+\/2 H 2+x/2 e -H/T for H > T > T*. These predic- 
tions seem to be in agreement with the experimental data 
in a series of different materials where disorder plays a 
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Figure f : Schematic phase diagram Txx where x controls the 
phase transition. Also shown is the behavior of the magnetic 
susceptibility, x(T), in the proximity of the quantum critical 
point. 

Using the Hertz theory of quantum critical phenom- 
ena || and treating the problem of an Ising field theory 
in the presence of a single impurity, Millis, Morr, and 
Schmalian (MMS) || have cast doubts on the possibility 
of Griffiths singularities in these materials based on the 
over-damping of the cluster (droplet) dynamics due to 
the electronic degrees of freedom. This effect is the re- 
sult of the slow decay of the droplet profile in the single 
impurity case. MMS have extended their calculation to 
a finite density of impurities 0] and conclude that over- 
damping destroys quantum Griffiths singularities leading 
to x{T) oc 1/T times ln(T) corrections at all T < ujq 
where loq is the high temperature cut-off of the theory 
above which the droplets are thermally activated. This 
result is in agreement with predictions in ref. that 
X(T) oc l/(Tln(l/T)) at low temperatures. MMS, how- 
ever, do not observe any crossover with T (as it was found 
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in ref. pj) and conclude that power law behavior should 
not be observed. 

In what follows we demonstrate the mathematical 
equivalence of the theory proposed in ref. and by MMS 
|2|. Both theories have a non-universal energy scale, 
T* = uj e~ C2 < uj , so that for T* < T < u one has 
x{T) oc T~ 1+x and therefore quantum Griffiths singular- 
ities. T* plays the role of an effective Kondo temperature 
of the droplet that is a non-universal quantity dependent 
on the constant Ci defined below. As we show below, 
C2 depends on the microscopic couplings of the material 
and changes from system to system. While MMS take 
C2 « 1 || leading to a narrow region of quantum Grif- 
fiths behavior, we have argued in ref. |l( that Ci ^> f 
leading to T* <C oj and therefore to a sizeable region in 
T where power law behavior should be observed. We be- 
lieve this result to be in agreement with experiments in 
various different materials^ @- In this paper we revisit 
the theory proposed by MMS and study the size of the 
Griffiths region (power law behavior) as a function of C2. 
We will show that this region shrinks very fast with the 
decrease of C2. 

The starting point is the Hertz action for a field theory 
with Ising symmetry S = So + Si, where: 
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here Eq is a characteristic energy scale (assumed to be 
of the order of the Kondo temperature of the system) , £ 
is a characteristic length scale (of the order of the lattice 
spacing), £ is the correlation length, c is a characteris- 
tic velocity, T is the damping coefficient, and cj)(r, r) is 
the order parameter (</>(k, a;„) its Fourier transform in 
momentum and Matsubara frequency). In principle the 
quantities that appear in the action have to be obtained 
from a microscopic theory or to be used as fitting param- 
eters to the experiments. 

The main difference between the two approaches dis- 
cussed in ref. and ref. Q is that ref. describes a soft 
spin model in the continuum while in ref. Q] a hard spin 
approach in a lattice is used. In the first case the ampli- 
tude of the droplets, <p , can fluctuate statistically while 
in the second case the amplitude of a cluster of N spins is 
fixed by the spin S of the spins in the cluster. Using the 
single impurity solution discussed in ref. Q as a varia- 
tional ansatz, MMS derived the probability distribution, 
P(R,(/)o) for a droplet of size R and amplitude <pQ, and 
the renormalized droplet tunneling splitting, Wtun- 

In order to make the connection between the differ- 
ent theories we notice that the number of spins in a 
droplet is approximately given by: / dr</>(r)/(47r£o/3) oc 
(/>o(i?/£o) 3 j where it is assumed that the droplet exists 
in a region of size R only. In what follows, instead of 



using MMS's reduced variables: y = R/£, / = 0o(£/£o)j 
and u = V 2 (£,/£,o) (Vq is the strength of the disorder) we 
define a new variable: 
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and new parameters: 
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where 



a(y) = 1 + 3/y 2 , 
c 7 - Eo/{6T), 
C 2 = l + ln[Efe*/(6c 2 )} 



(4) 



Notice that the parameters defined in J^J are indepen- 
dent of R and 4>o- Using these variables we can rewrite 
the main results of ref. [2j, namely, the droplet tunneling 
splitting (eq. (36) of ref. Q), 
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and the probability distribution {N[y 3 ,f 2 ] in eq. (15) of 
ref. 0), 



P(N,f)ocN 1 - e- N ^ N ^ 



where 6 — 5/2 and 



N 6 (f) = 



(6) 



(7) 



We can now make a direct comparison between the two 
different theories. Notice that the tunneling splitting is 
only dependent on the number of spins in the droplet, 
N, and it is identical to eq. (5.6) of ref. 0] if we assume 
a single cut-off energy scale (that is, set W — uo and 
(/3 = 1, see eq. (4.84) in ref. 0]). This result allows us 
to identify N c as the maximum number of spins in the 
droplet such that, above this number, tunneling ceases 
to occur, and v is the effective damping coefficient (the 
equivalent of 7 in eq. (5.6) of ref. [jj]). The probability 
distribution © is the equivalent of eq. (5.7) of ref. 0] ex- 
cept that in ref. |2j this probability distribution depends 
not only on the number of spins in the droplet but also 
on their amplitude. Although in a soft spin model there 
are fluctuations in the size of the droplet, those fluctua- 
tions are limited by the exponential term in ©. In fact, 
MMS find that the droplets that contribute most to the 
magnetic susceptibility are such that 
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If one replaces JSJ into eq. Q we find, for £ ^> £ , 
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In ref. it is shown that (see eq. (4.6) of ref. Q): 
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where D = 2.54 is the fractal dimension of the cluster in 
three dimensions. Therefore, © and (|10f> are very similar 
(the difference in the exponents is due to the different 
nature of the percolation in the two problems) and show 
that iVj diverges as £ — > oo. 

It should be clear at this point that the two main ele- 
ments of the theory, namely the tunneling splitting and 
its distribution, although calculated in very different way, 
are essentially the same in both theories. The main ques- 
tion is why MMS reach the conclusion that power law 
behavior, that is, quantum Griffiths singularities cannot 
be observed while in ref. Q it was found that there is 
a wide range in T where power law behavior should be 
observed. To answer this question let us consider the dis- 
tribution of tunneling splittings, that is, the probability 
of finding a droplet with tunneling splitting between A 
and A + dA, (instead of the distribution of cluster sizes). 
It is easy to convert from one to the other using J5J and 
© (we define A = cj t un): 



P(A) = P(N(A)) 
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It is straightforward to see that there is a characteristic 
tunneling splitting A c , 
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such that for A ^> A c one has: 

P(A) oc AVW- 1 



(12) 
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which leads to x{T) x T 1+A for T > A c in agreement 
with eq. (6.1) of ref. [l(. Using © and © we find: 
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which is eq. (44) of ref. 0. This last result can be com- 
pared with eq. (6.1) of ref. where: 
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and both theories predict that A — > as one approaches 
the QCP (see figDJ. This behavior is characteristic of 



quantum Griffiths singularities. Once again the differ- 
ence in the dependence of A with the dimensionality 
comes from the fact that ref. deals with a percola- 
tion problem in a lattice while ref. Q studies an impurity 
problem in the continuum. Nevertheless, the overall con- 
clusion that A vanishes at the critical point is common 
to both theories. On the other hand, if A <C A c we have 
from ifTTl) : 



P(A) oc 



Aln(cj /A) 
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which leads to the result that x(T) oc l/(TTn(l/T)) in 
agreement with eq. (5.18) in ref. |lj and with the con- 
clusion expressed in ref. |2j that the susceptibility is 1/T 
times logarithms. Thus, we have unequivocally shown 
that the theory presented by MMS has a crossover tem- 
perature T* = A c from power law to l/(Tln(l/T)) from 
high to low temperatures and alsoproved the equivalence 
of the results of ref. [l| and ref. [2| • 

The reason why the crossover is not discussed in ref. Q 
can be easily understood if we rewrite A c in l|12(l in terms 
of the parameters in 1(5}: 
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which means that the crossover depends on the value of 
C2, that is, depends on a non-universal number. By as- 
suming that c/£ » E MMS find C 2 » 1. This choice 
implies that A c /uo ~ 0.4 which is not particularly small 
and leads to a limited dynamical range for the power law 
to be observed (in Fig. 3 of ref.Q the plots stop at fre- 
quencies of order of O.lwo while the crossover should be 
observed at 0.4wo)- Thus, it is of great interest to under- 
stand how the magnetic response of the system depends 
on C2. 

We see from (@J that C2 depends on the details of the 
lattice through £0 and c as well as the Kondo energy 
scale through Eq. Any physical crossover, by definition, 
depends on non-universal parameters. A famous exam- 
ple is the Kondo problem where the Kondo temperature, 
Tk = ujoe^ 1 / 9 (where g is the strength of the spin ex- 
change interaction between localized moments and con- 
duction electrons), is also a non- universal function of the 
cut-off and, like the problem at hand, is an exponential 
function of effective coupling. As in the Kondo case the 
crossover temperature in this problem is exponentially 
sensitive on the choice of Ci. The crossover can be eas- 
ily visualized if one considers, for instance, the function 
J res (<x>) defined in ref. [2|. This function is directly related 
with the contribution to x(^) coming from the fluctuat- 
ing clusters through eq. (38) of ref. [2j: 



x(T) = £ 
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which is the analogue of eq. (5.15) of ref. 0. In fig[21 
we plot ac 2 Jres(w)/ Jres(0) versus uj/luq for Vq = 0.5, 
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Figure 2: Plot of J re s(w)/Jres(0) as a function of uj/ujo for 
C2 = 10 (open squares), C2 = 4 (closed triangles), C2 = 3 
(open circles), C2 = 2 (closed diamonds), C2 = 1 (open tri- 
angles), C2 = 0.1 (closed squares). The other parameters are 
given in the text. The crossover frequency from over-damping 
to under-damping is roughly the frequency associated with the 
maximum in each curve. 

c 7 =0.1 and £ = 20£o for various values of C2 where ac 2 
is a scaling number that allows the maxima in each curve 
to be at the same height («o.5 = 0.29, ai = 1, a.% = 3, 
as = 5.2, a4 = 7.4, aio — 15.5). One can clearly see 
that the crossover frequency, which is essentially marked 
by the maximum of each curve, occurs at low frequencies 
for C2 > 1 and at large frequencies when C2 < 1. This 
indicates that the region of power law behavior shrinks, 
that is, as the dissipation in the system increases, as 
C2 decreases. The dissipation depends strongly on the 
value of the coupling of the spins to the metallic envi- 
ronment, that is, on the Kondo exchange coupling Jk 
and therefore it depends on the microscopic details. 
These couplings change from system to system leading to 
different values of T*. It is very hard to determine the 
value of these couplings from first principle calculations 
and one has to rely on the experimental data in order 
to fit the value of C2 as one usually does in impurity 
systems in regards to the Kondo temperature. This re- 
sult explains why certain non- Fermi liquid materials show 
strong power law behavior than others 0, . 

We note that the conclusions of this work are rele- 
vant for systems with Ising symmetry. Recent results 
for particular realizations of Heisenberg quantum clus- 



ters [l(J and in studies of metallic Heisenberg magnets 
|l l| have shown that tunneling is not suppressed and that 
Griffiths-McCoy singularities are possible down to T = 0. 
Nevertheless, essentially all systems of experimental in- 
terest where non-Fermi liquid is observed have magnetic 
anisotropics due to spin-orbit coupling and crystal field 
effects |3J and therefore are in the Ising universality class 
discussed here. 

In summary, we have shown that the theory proposed 
by MMS in ref. 0] has generically the same crossover be- 
havior that we predicted in ref. where for T* < T < ojq 
quantum Griffiths singularities with non-universal expo- 
nents, x(T) oc T~ 1+A , should be observed and that for 
T < T* cluster freezing leads to x(T) oc l/(Tln(l/T)) 
(see fig'Q. The value of T* , however, varies from sys- 
tem to system and cannot be obtained within the theory. 
We have argued in ref. [l| T* <c an d therefore power 
law behavior should be clearly visible in a crossover re- 
gion. This conclusion seems to be in agreement with the 
experimental data in various materials 0, 0, 0] • 
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